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ABSTRACT 

We have identified galaxies near two quasars which are at the redshift of damped Lyman-a 
(DLA) systems in the UV spectra of the quasars. Both galaxies are actively forming stars. One 
galaxy has a luminosity close to the break in the local galaxy luminosity function, L* , the other is 
significantly fainter than L* and appears to be interacting with a nearby companion. Despite the 
strong selection effects favoring spectroscopic identification of the most luminous DLA galaxies, 
many of the spectroscopically- identified DLA galaxies in the literature are sub-L*, suggesting 
that the majority of the DLA population is probably sub-L*, in contrast to Mgii absorbers at 
similar redshifts whose mean luminosity is close to L* . 

Subject headings: Quasars: absorption lines — galaxies: distances and redshifts — quasars: individual 
(FBQS J0051+0041, FBQS J1137+3907) 



1. Introduction 

Damped Lya (DLA) absorbers in quasar spec- 
tra are the class of absorbers with the highest 
column density of neutral hydrogen {Nhi > 2 x 
lO^^cm^^). They are important as it can be shown 
that they contain the bulk of the neutral hydrogen 
content of the Universe (e.g. Storrie-Lombardi & 
Wolfe 2000). The similarity of the column densi- 
ties of DLAs to those through the disks of spiral 
galaxies have lead some authors to consider them 
as the direct progenitors of present-day spirals 
(e.g. Prochaska & Wolfe 1998). However, dwarf 
galaxies and low surface brightness galaxies (LS- 
BGs) also can have similar columns (Boissier, Per- 
oux & Pettini 2003), and metallicity arguments 
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support dwarf galaxies (Pettini et al. 1999). The- 
orists have interpreted them in the context of hi- 
erarchical galaxy formation models as sub-units 
which will eventually merge to form L* galaxies 
today (Haehnelt, Steinmetz & Ranch 1998; Mailer 
et al. 2001) or as dwarf galaxies in which star for- 
mation is supressed by supernova feedback (Efs- 
tathiou 2000). 

The reason for the proliferation of theories for 
the nature of DLAs is that reliably identifying the 
galaxies responsible for the absorption, even at 
z < 1, has proven to be very diffictdt. Imaging 
studies of z < 1 absorbers, e.g. Lc Brun et al. 
(1997); Rao et al. (2002) indicate that a range 
of galaxy types might be responsible for produc- 
ing DLAs, but few candidate DLAs have mea- 
sured redshifts, so reliable constraints on the na- 
ture of the DLA population remain poor. Para- 
doxically, Mgii absorbers, with much lower col- 
umn densities on average, have been easier to iden- 
tify. They seem to be sampling the normal galaxy 
population, with a mean luminosity of ^ L* at 
< z >= 0.6 (Steidel, Dickinson & Pcrsson 1994; 
Steidel et al. 2002). Possible reasons for the diffi- 
culty of identifying the galaxies causing DLA sys- 
tems include small impact parameters, low stellar 



1 



luminosities or low surface brightnesses. So far 
only six 2 < 1 DLAs have been spectroscopically- 
confirmed (and all of these are a.t z < 0.6). Several 
arc dwarf or low surface brightness galaxies, the 
best-studied example of which is SBS 1543+5921 
(Bowen, Tripp & Jenkins 2001), but the DLA ab- 
sorber of AO 0827+24 is a fairly normal, lumi- 
nous spiral (Steidel et al. 2002). Many more sug- 
gested identifications have been made on the basis 
of proximity to the quasar and consistent photom- 
etry, but the reliability of these is unclear. 

DLAs are easy to identify in the spectra of high 
redshift {z ^ 2) quasars as they show up clearly 
in observed optical spectra, but they are hard to 
identify at lower redshift. By adopting the strat- 
egy of taking UV spectra with the Hubble Space 
Telescope {HST) of the highest rest-frame equiv- 
alent width {EWo) Mgii absorbers in the litera- 
ture, Rao & Turnshek (2001) were able to identify 
a sample of lower redshift {z ~ 0.3—1.7) DLAs (see 
also Boisse & Bergeron 1998). High EWq Mgii 
absorbers also typically show high EWo Fell ab- 
sorption lines, and these have been hypothesised 
to be better tracers of cool gas than Mgii (Berg- 
eron & Stasihska 1986). We have therefore used 
the FIRST Bright Quasar Survey (FBQS; White 
et al. 2000; Becker et al. 2001) to select a simi- 
lar sample with EWq (Fell 2600A)> lA. We have 
taken spectra of 26 of the sample with STIS on 
the HST to look for DLA systems. We have fol- 
lowed up this sample using the Keck telescope, 
taking near-infrared images in sub-arcsecond see- 
ing, medium-resolution spectra of the quasars, 
and low-resolution spectra of the candidate ab- 
sorbers. We are mid-way through spectroscopic 
follow-up of the FBQS sample, and have iden- 
tified four galaxies with redshifts the same as 
that of the Fell absorbers. Of these four, how- 
ever, only two are associated with absorbers which 
are definitely damped. These two galaxies, near 
the quasars FBQS 0051+0041 (R.A. 00'^51'"30!49 
Dec. +00°41'"49'.'91 (J2000); absorber redshift, 
Zabs = 0.74; quasar redshift, zq = 1.19; appar- 
ent magnitude B = 18.5) and FBQS 1137+3907 
(R.A. ll'^37™09!46 Dec. +39°07™23'.'60 (J2000); 
Zabs = 0.72; ZQ = 1.02; B = 18.4), form the sub- 
ject of this paper. 

We assume a cosmology with = 0.3, a = 
0.7 and Ho = 70kms"^Mpc"^ throughout. Where 
we have referenced galaxy luminosities, L, to the 



present-day values of the break in the luminosity 
function, L* we have used L* values of Norberg et 
al. (2002) in B, Lin et al. (1996) in R and Loveday 
(2000) in K, adjusted to our passbands and cos- 
mology as necessary (M^ = -20.2, = -21.1, 
= —24.4). The fc-corrections have been de- 
rived using a model 5Gyr old Sd galaxy from the 
Pegase library (Rocca-Volmerange & Fioc 1997). 
We picked this model as representative of a mod- 
erately young, star-forming galaxy, a decision mo- 
tivated by the discovery of [On] emission from the 
candidate DLA galaxies. 

2. Imaging Observations 

FBQS 0051+0041 was imaged with the Echel- 
lette Spectrograph and Imager (ESI) on Keck- 
II and FBQS 1137+3907 with the Near-Infrared 
Camera (NIRC) on Keck-I. Details of the observa- 
tions are given in in Table 1. A sequence of short 
(3-5min) dithered exposures of each object were 
taken. After dark subtraction and flat-fielding 
with dome flats the data were analysed using dim- 
sum and XDIMSUM in IRAF, respectively. The 
DIMSUM package performs two passes on the data. 
In the first pass, a background image is made from 
a median of the dithered frames and the data are 
combined with integer pixel shifts to make masks 
which are used to mask out objects when con- 
structing an improved sky background image for 
the second pass. In the second pass, the improved 
background is subtracted from the images before 
they are shifted and coadded, and a magnification 
factor may be applied to the data. The pixel scale 
of the ESI image was suSiciently small compared 
to the seeing that no magnification was used in the 
second pass for FBQS 0051+0041. The operation 
of the DIMSUM package is explained in more detail 
in Stanford, Eisenhardt & Dickinson (1995). The 
XDIMSUM package is a modification of dimsum by 
L. Davis at NO AO. Amongst other improvements 
it runs significantly faster and allows non-integer 
magnification factors in the second pass. We used 
a factor of 2.3 for the FBQS 1137+3907 image. 

Photometry for FBQS 0051+0041 was obtained 
using observations of the Landolt (1992) standard 
107-215, and that for FBQS 1137+3907 from ob- 
servations of the standard star P266-C from Pers- 
son et al. (1998). Both quasar fields were imaged 
in photometric conditions. The 5ct limiting magni- 
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Fig. 1.— The field of FBQS 0051+0041 in /. The 
quasar has been partially subtracted from near the 
center of the image, and its position is marked 
with a cross. The DLA candidate is galaxy 'a' 
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Fig. 2.— The field of FBQS 1137+3907 in K' . 
The quasar has been partially subtracted from 
near the center of the image, and its position is 
marked with a cross. Galaxies 'a' and 'c' are most 
probably associated with the DLA system. 



tudes in 2-arcsec diameter apertures are / = 25.4 
for the FBQS 0051+0041 field and K' = 21.9 for 
the FBQS 1137+3907 field. 

We performed approximate subtractions of 
the quasars using point spread functions de- 
rived from either from a star in the same im- 
age (FBQS 0051+0041) or from other quasars 
observed on the same night (FBQS 1137+3907). 
The PSF was fitted with a least-squares fit us- 
ing code developed for fitting quasar host galaxies 
(Lacy et al. 2002). The particular algorithm em- 
ployed here subtracted the PSF such that the 
mean value within a radius rin from the center 
of the quasar was set to zero. The values of nn 
used were O'.'S for FBQS 0051+0041 and 0'.'3 for 
FBQS 1137+3907. Neither PSF subtraction was 
very satisfactory, as we made no serious attempt 
to obtain good, unsaturated PSFs for the quasars. 
They have successfully removed the wings of the 
PSFs, however. The resulting images are shown 
in Figures 1 and 2. 

3. HST Spectroscopy 

The quasars were observed with the STIS with 
the MAMA detectors and the G230L grating as 
part of program 9051. The 0'.'2 slit was used, 
giving a spectral resolving power of 330 at the 
redshifted wavelength of Lya, corresponding to a 
velocity resolution of f« 900kms~^. Each quasar 
was observed for two 1100s exposures, the sec- 
ond exposure being dithered 2 along the slit. 
The UV spectrum of FBQS 1137+3907 shows 
it to be a high-ionization broad- absorption line 
(BAL) quasar, with a BAL trough blueshifted by 
13000kms~^. Fortunately, the Lyman-a absorp- 
tion line is clear of the BAL troughs, falling just 
to the red of the Ovi trough in the quasar spec- 
trum. The spectra are shown in Figures 3 and 
4. 

A detailed discussion of the HST spectra and 
measurements of the damped absorbers are pro- 
vided in a companion paper (Storrie-Lombardi, 
et. al. 2003, in preparation). The column densities 
of the damped absorbers were initially estimated 
from measurements of the equivalent widths of the 
absorption lines with the redshift determined from 
the low ionization metal lines from our ground- 
based spectra. Profiles were then fit iteratively to 
determine the best fit column density and redshift 
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for HI absorption feature. The column densities 
should be accurate to ±0.1 and the redshifts to 
±0.005. Though higher resolution spectra cover- 
ing multiple lines in the Lyman series would be 
ideal for this purpose, the substantial HST time 
that would be involved in obtaining them makes 
such observations impractical for large samples of 
objects. At higher redshifts, we have shown via ob- 
servation and simulation that moderate resolution 
spectra provide an excellent estimate for the mea- 
sured HI column densities (e.g. Storrie-Lombardi 
et. al. 1996; Storrie-Lombardi & Wolfe 2000), 
which gives \is high confidence that the features 
in these spectra are damped and are not blends of 
lower column density lines. 

4. Keck spectroscopy 

The candidate absorbers were observed with 
the Low Resolution Imaging Spectrometer (LRIS) 
on Keck-I using the 300 line mm^-*^, 5000 A blaze 
grating in the red arm and a 1-arcsecond slit (see 
Table 2). This combination results in a spectral 
resolving power of « 500, corresponding to a ve- 
locity resolution of « 600kms~^. The wavelength 
range covered was from w 5000 - lOOOOA. 

In addition, observations of the quasars were 
made with ESI on Keck-II to obtain medium- 
resolution spectra of the absorption lines. The Fell 
2586 and 2600A lines fell in a noisy region near an 
order overlap in the ESI spectra, but the Mgll 2796 
and 2803A lines were detected at high signal-to- 
noise in both quasars. FBQS 0051+0041 was ob- 
served with ESI on 2000 December 31 with a 0'.'75 
slit, and FBQS 1137+3907 on 2001 March 20 with 
a 1 slit, yielding resolving powers of 5400 and 
4100 respectively. The ESI spectra of the quasars 
in the region around Mgii are shown in Figures 5 
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Fig. 4. The STIS spectrum of FBQS 
1137+3907. The upper plot has the quasar Lya 
emission (at zq = 1.02) and BAL features labelled, 
the lower plot has absorption features associated 
with the DLA system (zabs = 0.72) labelled. 
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Fig. 5. The Keck ESI spectrum of FBQS 
0051+0041 in the region of the Mgii absorption 
lines. 
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Fig. 3.— STIS spectrum of FBQS 0051+0041. 
The quasar {zq = 1.19) and absorber (^abs = 0.74) 
features are labelled. 
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Fig. 6.— The Keck ESI spectrum of FBQS 
1137+3907 in the region of the Mgii absorption 
lines. 
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and 6, and the LRIS spectra of the galaxies at the 
absorber redshifts in Figures 7 and 8. 

5. The absorber of FBQS 0051+0041 

The HST spectrum shows that this absorber 
is a borderline DLA, with a column density just 
above the division between Lyman-limit and DLA 
systems (Table 3). Similarly, the EWo of the Fell 
absorption only just satisfies our selection criteria. 

The candidate absorber galaxy, galaxy 'a' in 
Figure 1, seems, at first sight, to be a fairly nor- 
mal, star forming galaxy, offset by 24kpc from the 
quasar (Figures 1 and 7). The spectrum, however, 
shows a small mis-match between the [On] 3727 
redshift of 0.740 and the absorption line rcdshift 
from the Ca II H+K lines (0.738), though the 
signal-to-noise ratio of the absorption lines is low. 
Two distinct velocity components are also seen in 
the low-ionization Mgll, Mgl and Fell quasar ab- 
sorption lines in our ESI spectrum, one a.t z = 
0.7393 and one at z = 0.7402 (i.e. a IGOkms"! ve- 
locity difference). Whether these two components 
correspond to two galaxies, the candidate and a 
dwarf companion nearer the quasar, or reflect in- 
ternal gas motions in the candidate absorber is 
unclear. 

6. The absorber of FBQS 1137+3907 

This is a much higher column density absorber 
than that of FBQS 0051+0041. The HST spec- 
trum shows several strong UV metal absorption 
lines in addition to the damped Lya line (Figure 
4). 

The LRIS slit was placed to cover both candi- 
date absorbing galaxies, a and c, and line emission 
is seen from both of them (barely resolved from 
each other) at the wavelength corresponding to 
[Oil] 3727 at z = 0.7185+0.0003 (Figure 9). There 
is also a marginal detection of H/3 at the same red- 
shift (Figure 8). Although we sec no other features 
from a or c, it seems unlikely that the strong emis- 
sion line can be anything other than [On] at the 
absorber redshift, given its high equivalent width 
and the lack of lines nearby in the spectrum. We 
also took a spectrum with the slit aligned with 
object 'b'. This showed that 'b' had blue optical 
continuum emission, but no strong features from 
which to obtain a redshift. As it is further from 
the quasar than 'c' and has no features at wave- 



lengths corresponding to the absorption rcdshift 
we assume that it is unrelated to the DLA sys- 
tem. 
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Table 1: Imaging of the FBQS DLA fields 



Quasar field 


Telescope/Instrument 


Filter 


Observation 


Total integration 


Seeing 








date (UT) 


time (minutes) 


(arcsec) 


FBQS 0051+0041 


Kcck-II/ESI 


/ 


2001 Juno 22 


20 


0.8 


FBQS 1137+3907 


Keck-I/NIRC 


k' 


2002 May 25 


27 


0.4 



Table 2: Observed properties and spectroscopic observations of the FBQS DLA galaxies 



Quasar field 


^gal 




Mag. 


Impact 

parameter 
(arcsec) 


Impact 

parameter 
(kpc) 


Abs. 
mag. 


AMfc 


L/L' 


UT Dale 
(dcg) 


Slit PA 


Exposure 
time (s) 


FBQS 0051+0041a 


0.740 


/ 


= 22.0 


3.3 


24 


Mb = -20.0 


-0.11 


0.8 


2001 Sopt 25 


26 


2 X 1200 


FBQS 0051+0041b 


7 


/ 


= 23.3 


4.6 


7 














FBQS 0051+0041C 


0.070 


/ 


= 23.2 


7.2 


9.6 


Mr = -13.9 


-0.15 


0.001 


2001 Sept 25 


26 


2 X 1200 


FBQS 1137+3907a 


0.7185 


K 


= 19.5 


2.5 


18 


Mk = -23.0 


-0.08 


0.3 


2002 Dec 31 


126 


2 X 1800 


FBQS 1137+3907b 


7 


K 


= 20.4 


1.8 


7 








2002 Dec 31 


228 


2 X 1800 


FBQS 1137+3907C 


0.7185 


k' 


' = 19.8 


1.5 


11 


Mk = -22.8 


-0.08 


0.2 


2002 Doc 31 


126 


2 X 1800 



Notes: AM^ is the amount that was added to the absolute magnitudes to account for the difference in the 



flux density of the galaxy between the observed band and the band the absolute magnitude was calculated 
in. This correction depends only on the assumed galaxy SED and reduces to the usual ^-correction when 
the two bands are the same. 



Table 3: Details of the DLA systems 



Quasar 


^abs 


Mgll EWo 


Fell EWo 


\g{Nm 






(2796) 


(2600) 


(cm2)) 


FBQS 0051+0041 


0.740 


2.4 


1.0 


20.4 ±0.1 


FBQS 1137+3907 


0.719 


3.0 


2.5 


21.1 ±0.1 



Tabic 4: Physical properties of galaxies with spectroscopic rcdshifts consistent with Zabs 



Absorber 


^gal 


lg(NHI 

(cm=)) 


Impact 


Observed 


Absolute 


AMfc 


L/L' 


Type 


Refs. 


candidate 




parameter 
(kpc) 


mag. 


mag. 










SBS 1543+5921 


0.009 


20.4 


0.4 


R = 16.3 


Mr = -16.7 




0.02 


LSBG 


9,10 


OI 363 (B0738-I-313) Gl 


0.221 


20.9 


20 


R = 20.8 


Mr = -19.4 


0.2 


0.2 


dE 


3,4 


FBQS 1137+3907C 


0.720 


21.1 


11 


k' = 19.8 


Mk = -22.8 


-0.1 


0.2 


spiral? 


1 


PKS 1127-145 


0.312 


21.7 


<7 


R = 21.6" 


Mb = -19.3 


-0.2 


0.4 


LSBG 


7,8 


PKS 1243-072A 


0.437 


fa 20.5'' 


12 


I = 21.3 


Mr = -20.1 


0.2 


0.4 


starburst/Seyfert 2 


8 


FBQS 0051-l-0041a 


0.740 


20.4 


24 


I = 22.0 


Mb = -20.0 


-0.11 


0.8 


spiral? 


1 


B2 08274-243 Gl 


0.526 


20.3 


36 


/ = 20.3 


Mr = -21.6 


0.2 


1.6 


spiral 


5,6 


AO 0235-M64 A" 


0.524 


21.7'' 


13 


m785 > 20.3 


Mr > -22.6 


0.3 


< 4 


Seyfert 1 


2 



Notes: " object Al in ref. 2 is closer to the quasar but has no spectrum, if it is at the absorber redshift it has 
{L/L*)b = 1-0, A has a bright nuclear component, and the luminosity of the underlying galaxy is uncertain; 
we have quoted upper limits based on the sum of the host and AGN fluxes. column density from 21cm 
line absorption, dependent on the assumed spin temperature. ^ sum of components 1-4 of Rao et al. (2003). 
See Table 3 for the definition of AMfc. References: (1) this paper; (2) Burbidge et al. 1996; (3) Turnshek et 
al. 2001; (4) Cohen 2001; (5) Steidel et al. 2002; (6) Rao et al. 2003; (7) Lane et al. 1998; (8) Kanekar et al. 
2002; (9) Reimers & Hagen 1998; (10) Bowen et al. 2001 
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The ESI spectrum shows a comphcatcd veloc- 
ity profile, with two to four distinct low-ionization 
absorption components seen, depending on which 
line is examined. The two deepest features, visible 
in all the lines, are at z = 0.7184 and z = 0.7193 
(i.e. a 160kms~^ velocity difference). Thus the 
[On] emission is most likely to be associated with 
the deeper, less redshifted absorption component. 

7. Discussion 

The DLA absorber candidate close to FBQS 
0051+0041 has a luminosity of approximately L* , 
whereas that associated with FBQS 1137+3907 is 
about 0.2L*. The galaxy types are unclear. Both 
are apparently forming stars, possessing both blue 
continua and [On] emission lines, and the com- 
plicated velocity structures of the Mgil absorp- 
tion lines seen in the ESI spectra, might be in- 
dicative of starburst-induced winds (e.g. Ellison 
et al. 2003). Further evidence for this is the offset 
between the emission and absorption redshifts in 
FBQS 0051+0041a, if real. Thus there is some ev- 
idence that both may be starburst galaxies. How- 
ever, the continuum slope of FBQS 0051+0041a 
is a little too red unless the starburst is reddened 
by dust, and FBQS 1137+3907c has a high axial 
ratio, suggesting a spiral rather than an irregular 
galaxy (though the shape of 'c' may also be due to 
it being composed of a chain of knots). Higher spa- 
tial resolution multicolor imaging will be required 
to determine the nature of these galaxies. 

Table 4 summarizes the properties of the eight 
known candidate DLA galaxies at z < 1 which 
spectroscopy has shown to lie at the absorber red- 
shift. Selection effects will mean that the absorber 
population for which spectroscopic identification 
is possible will be biased. Observational require- 
ments for successful spectroscopy, namely an im- 
pact parameter ^ 1 arcsec and an apparent mag- 
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Fig. 7. — The spectrum of galaxy 'a' {z = 0.740) 
near FBQS 0051+0041. 
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Fig. 8. — The sum of the extracted spectra of 
galaxies 'a' and 'c' near FBQS 1137+3907. 
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Fig. 9. — The 2D spectrum of galaxies 'a' and 'c' 
near FBQS 1137+3907 in the region of the [On] 
emission line in the galaxies. The trace of the 
quasar was averaged for 30 pixels either side of 
the region shown and the mean subtracted from 
this portion of the image. The center of the un- 
subtracted quasar trace is indicated by the dashed 
Une at arcsec. The seeing FWHM was 0'.'92 (4.3 
pixels), so 'a' and 'c' are only marginally resolved. 
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nitudo brighter than R ^ 24, wiU naturaUy fa- 
vor massive galaxies. For our FBQS sample dis- 
cussed in the introduction, and from which FBQS 
0051+0041 and FBQS 1137+3907 were drawn, 
we have so far only successfully identified candi- 
dates and obtained spectra with redshifts consis- 
tent with those of the absorber system for four out 
of thirteen of our high EWo absorption systems. 
In addition, there is another subtle selection ef- 
fect from the way most (though not all) of the 
DLAs in Table 4 have been selected, namely on 
the basis of high EWq metal absorption lines. This 
will tend to bias the sample towards higher mctal- 
licity, more massive, galaxies. Working against 
these two effects is the possibility that many of 
the more massive DLA galaxies arc dusty enough 
to significantly extinguish the light from the back- 
ground quasar, and cause the quasar to drop out 
of optically-selected samples. This may tend to 
produce a bias towards relatively dust-free, less 
massive hosts, and may explain why the Mgll ab- 
sorbers, with their generally higher impact pa- 
rameters and lower Hi column densities, are ob- 
served to have luminosities around L* . However, 
the results of the survey of radio-loud quasars of 
Ellison et al. (2001) suggest that this is not an 
important effect. Any such effect will also be 
much reduced for FBQS quasars as they arc se- 
lected on the Palomar i^-plates with a relatively 
red {O — Er^B — R< 2.0) color selection cri- 
terion. Thus, despite most of the selection effects 
favouring successful spectroscopic identification of 
the most luminous absorbing DLA galaxies, it is 
clear from Table 4 that many are sub-L*, and this 
suggests that the majority of the DLA population 
is probably sub-L*. 

Both our DLA candidates appear to be ac- 
tively forming stars, consistent with them being 
relatively gas-rich galaxies. A large fraction of 
the known spectroscopically-confirmed DLA sys- 
tems at low redshift listed in Table 4 appear to be 
merging or interacting in some way, for example, 
FBQS 1137+3907C, PKS 1127-145 and PKS 1243- 
072 A (though obviously the statistics are poor at 
present, and few apparent pairs have been spec- 
troscopically confirmed). This is intriguing in the 
light of theories in which DLA galaxies are dwarf 
galaxies which merge to become massive galaxies 
today (Haehnclt et al. 1998; Mailer ct al. 2001). 
Though these theories are more applicable to the 



high redshift systems, the lower redshift DLAs 
may represent the tail end of this activity. 
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